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На основе реологической модели В.Н. Покровского, описывающей поведение наполненных полимерных систем, рассматрива- 
ется задача о стекании пленки нелинейной вязкоупругой жидкости по боковой поверхности вертикального цилиндра. Приво- 
дится математическая постановка задачи, методика ее решения, обсуждаются результаты расчета и анализируются факторы, 
оказывающие влияние на процесс стекания пленки. 


Введение 

Пленочные течения реологически сложных 
жидкостей встречаются во многих отраслях совре- 
менного производства - химической технологии, 
энергетике, транспортировке нефтепродуктов и др. 
Для нанесения на детали и изделия различных по- 
крытий и смазок используют методы макания и 
окунания. При быстром опорожнении сосудов и 
резервуаров различного назначения на их стенках 
остаются пленки жидкости, постепенно стекаю- 
щие вниз. Поэтому важно определить характери- 
стики таких пленок (толщину, расход, время стека- 
ния, оставшееся на стенке покрытие) в зависимо- 
сти от реологических характеристик жидкости [1]. 

В большинстве случаев реальные материалы, 
находящиеся в вязкотекучем состоянии, обладают 
ярко выраженными неньютоновскими свойствами, 
такими, как нелинейная вязкость, пластичность и 
вязкоупругость. Неправильный учет влияния рео- 
логических факторов на движение пленок при эк- 
сплуатации различных технологических устройств 
приводит к существенным ошибкам в расчете. 

При теоретическом описании технологических 
процессов переработки наполненных полимеров 
встает вопрос о локальных законах течения, учиты- 
вающих существенные реологические особенности 
материала. Имеющиеся теории течения часто осно- 
вываются на модели нелинейной вязкой или вязко- 
пластичной жидкости. Такое описание часто оказы- 
вается недостаточно полным и не включает очень 
существенных, обнаруженных на опыте, нелиней- 
ных эффектов, связанных с наличием релаксацион- 
ных процессов в жидкости, таких, как появление 
нормальных напряжений при сдвиге и возникнове- 
ние остаточных напряжений при остановке тече- 
ния, разбухание струи при выходе из формующей 
головки и навивки жидкости на вращающейся 
стержень (эффект Вейсенберга) и др. [2-5]. 

1. Постановка задачи 

Для описания поведения реологически сложных 
жидкостей не существует универсальных определяю- 
щих уравнений, пригодных, с одной стороны, для 
всех нелинейных материалов, и с другой - удовлетво- 
рительно работающих в широком диапазоне скоро- 
стей сдвига. Нелинейные эффекты определяются спе- 


цификой материала, и это надо учитывать при форму- 
лировке реологических уравнений динамики жидко- 
сти. Одной из реологических моделей, достаточно хо- 
рошо описывающей поведение наполненных поли- 
мерных систем на основе линейных полимеров, явля- 
ется модель, предложенная В.Н. Покровским [6]. 

В соответствии с этой моделью система уравне- 
ний движения линейных полимеров, кроме обыч- 
ных уравнений движения в напряжениях и уравне- 
ния неразрывности, содержит совокупность ура- 
внений, состоящую из определения девиатора тен- 
зора напряжений 
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где X/ - декартовы координаты, ^ - внутренний 
параметр, определяющий меру деформации макро- 
молекулярных клубков, і - время, ѵ, - скорость, 
8 у - единичный тензор, <у і} - девиатор тензора на- 
пряжений, к - текущее значение индекса, по кото- 
рому выполняется суммирование, г и ц- время ре- 
лаксации и сдвиговая вязкость, которые предпола- 
гаются зависящими от приложенных напряжений, 
причем выполняется соотношение [6] 

Л-=Л = /т О = о и + Ър, ( 3 ) 

Ло т о 


где ц и г/о - начальные вязкость и время релаксации 
при /)=0, р - давление. Реологические свойства жид- 
кости в этой модели определяются значениями пара- 
метров г 0 и г/о, а также некоторой монотонно убываю- 
щей функцией//)), которая определяется по экспери- 
ментальным данным и для линейных полимеров мо- 
жет быть аппроксимирована, например, функцией [6] 

/(/)) = Ѵ(1 + Щ\ (4) 

где к - коэффициент, изменяющийся в интервале 
0,1. ..0, 2; п - показатель нелинейности. 

Рассмотрим задачу о стекании пленки нелиней- 
ной вязкоупругой жидкости, удовлетворяющей мо- 
дели В.Н. Покровского, по боковой поверхности 
вертикально расположенного прямого кругового 
цилиндра (рис. 1). 
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определяющих уравнений, полученных на основе 
модели В.Н. Покровского, ур. (1—4), и связываю- 
щих компоненты девиатора тензора напряжений с 
градиентом скорости, которые в данном случае 
принимают вид: 

^ 2Г =0о/ (1 + ^оА7оО"-^> ( 7 ) 
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Рис. 1. Г еометрическая схема течения пленки: г , г ~ цилин- 
дрические координаты; /? с - радиус цилиндра; 8 ~ 
толщина стекающего слоя жидкости; К § =К с +8 ~ рас- 
стояние до свободной поверхности 

Полагаем, что течение пленки является устано- 
вившимся, ламинарным, и безволновым. Боковая 
поверхность цилиндра поддерживается при по- 
стоянной температуре /=/ 0 , а температура окружаю- 
щей среды і=і ж . В этом случае необходимо учиты- 
вать кривизну орошаемой стенки. Ясно, что капил- 
лярные силы, обусловленные поперечной кривиз- 
ной, малы для цилиндров большого радиуса. Ван 
Россум [7] установил, что если число Гоугера, кото- 
рое определяется как 

у 

Со = К(р 8 /(2ст)У (5) 


где а - коэффициент температуропроводности 
жидкости. 

Уравнения (6-9) решаются при условии, что на 
боковой поверхности цилиндра выполняются 
условия прилипания, а температура стенки по- 
стоянна 


Ѵ 7 ~ О 

і = і 0 


при г = Я с , 



а на свободной поверхности пленки (г=К 8 ) задаются 
условия отсутствия касательного напряжения <х, и 
теплового потока д, определяемого по закону Фурье 


сг = О, 
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больше 1,8, то капиллярными силами в расчетах 
можно пренебречь. 

В соотношении (5) р - плотность жидкости; § - 
ускорение силы тяжести; сг - коэффициент по- 
верхностного натяжения. 

Из общих соображений следует, что при боль- 
ших значениях обобщенного числа Прандтля, ха- 
рактерных для всех реологически сложных жидко- 
стей, вблизи орошаемой поверхности должен обра- 
зовываться эффективный температурный слой, где 
происходит основное изменение температурного 
поля. В этом слое температура вдоль пленки изме- 
няется гораздо медленнее, чем в ее поперечном на- 
правлении. Поэтому при решении системы уравне- 
ний, описывающих процесс отекания пленки в 
условиях теплообмена, воспользуемся гипотезой 
квазистационарности, т. е. будем считать, что каж- 
дому мгновенному распределению температур в 
данный момент времени соответствует свое ста- 
ционарное распределение скоростей. Такое допу- 
щение позволяет решать гидродинамическую и те- 
пловую задачу раздельно. 

Таким образом, с учетом сделанных предполо- 
жений, для решения поставленной задачи необхо- 
димо решить систему уравнений, записанных в ци- 
линдрической системе координат и состоя- 

щую из уравнения движения 

г'^ (г О + РЯ = 0, (6) 


Выбирая в качестве масштаба длины радиус ци- 
линдра - г=К с , скорости - среднерасходную ско- 
рость стекания — ѵ= II и температуры - разность 
между температурой стенки цилиндра / 0 и началь- 
ной температурой жидкости і ж , запишем исходную 
систему уравнений (7-9) и граничные условия 


(10, 11) в безразмерной форме: 

-У(гаД + Ке/Рг = 0, ( 12 ) 
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ѵ 2 =0, 0 = 1 при г = 1, (16) 

сг_ г = О, =0приг = 1 + Д, (17) 

дг 


где Ѳ = — _ безразмерная температура; А =8/К с - 

^0 _ 1 ж 

безразмерная толщина пленки; Ке=рШ с /гі {) - число 
Рейнольдса; Рг= Л 2 /(§Я С ) - число Фруда; УІо=цЛ/К с - 
число Вейсенберга; Р о=Ш с /а - число Пекле. 
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Отношение числа Рейнольдса к числу Фруда, 
входящее в (12), представляет важный с гидродина- 
мической точки зрения критерий подобия, кото- 
рый обозначается через ^=Ке/Рг [8]. 

Найдем решение гидродинамической задачи. 

Нетрудно показать, что, используя уравнения 
(12-14) и первое из граничных условий (17), можно 
найти выражение для касательного напряжения <х,. 

а. г =^[(1 + Д) 2 -г 2 ] (18) 

2 г 

и получить связь между градиентом скорости и 
компонентой тензора напряжений а 1Г в виде 

7^=^(1 + 2ШеѴ)". (19) 

ОТ 



Рис. 2. Зависимость толщины слоя от числа И/е при различ- 
ных значениях И/; 1) 5; 2) 50; 3) 100 


Подставляя выражение (18) в (19), получим не- 
линейное дифференциальное уравнение первого 
порядка для определения профиля скорости по 
толщине пленки. Поскольку является функцией 
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В виду нелинейности дифференциального ура- 
внения (20) получить его аналитическое решение 
не представляется возможным. Поэтому это ура- 
внение решается численно с использованием пер- 
вого из граничных условий (16). Уравнение энер- 
гии, записанное в виде (15), также решается чи- 
сленно с использованием поля скорости, получен- 
ного из решения уравнения (20). 


2. Результаты расчета 

Проведенные численные эксперименты позво- 
лили выяснить влияние параметров течения и рео- 
логических характеристик жидкости на профили 
скоростей и температур, а также установить их 
влияние на толщину стекающей пленки А при за- 
данных значениях объемного расхода. Для опреде- 
ления величины А строился итерационный про- 
цесс с использованием соотношения 

і+д 

0^ = 2п \ гѵ 2 (г) сіг. 


На рис. 2 представлена зависимость толщины 
стекающего слоя вязкоупругой жидкости от числа 
Вейсенберга при различных значениях параметра \У. 

Из представленного рисунка видно, что толщи- 
на пленки уменьшается как с увеличением числа 
ЛЛ^е, так и с увеличением параметра \У. При этом с 
ростом \Уе толщина пленки стремится к некоторо- 
му стационарному значению. 


Зависимость значения максимальной скорости 
в том же диапазоне изменения параметров ѴеиѴ 
показана на рис. 3. 


^ шах 



Рис. 3. Зависимость максимальной скорости от числа И/е 
при различных значениях И/: 1) 5; 2) 50; 3) 100 

Представленные результаты расчета показыва- 
ют, что если толщина пленки с ростом параметров 
ДУе и ^ падает, то максимальное значение скоро- 
сти, которое достигается на свободной поверхно- 
сти пленки, растет. Рис. 4 иллюстрирует распреде- 
ление скорости по толщине пленки в зависимости 
от числа Вейсенберга (\Уе). 


ѵ 



Рис. 4. Распределение скорости по толщине пленки в зави- 
симости от И/е; 1) 50; 2) 100; 3) 150; 4) 200 
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Как видно из рисунка, повышение числа Вейсен- 
берга приводит к заметному увеличению скорости. В 
то же время представленный рисунок еще раз пока- 
зывает, что с ростом \Уе толщина стекающего слоя 
уменьшается. Распределение скорости для вязко- 
упругой жидкости зависит и от геометрических раз- 
меров или радиуса цилиндра, по которому стекает 
жидкость. Как показали численные эксперименты, 
увеличение радиуса цилиндра при условии сохране- 
ния постоянного значения объемного расхода приво- 
дит к понижению максимальной скорости и к умень- 
шению толщины стекающей пленки. 

В отличие от нелинейной вязкой и вязкопластич- 
ной жидкости при течении вязкоупругой жидкости 
нормальная компонента девиатора тензора напря- 
жений не равна нулю. При этом величина напряже- 
ния существенным образом зависит от числа \Уе: 

сг__ = 2ДѴе сг 2 . 

22 2Г 

На рис. 5 показаны результаты расчета нор- 
мального напряжения <х. по толщине пленки в за- 
висимости от числа \Ѵе. Из представленного рисун- 
ка следует, что распределение нормального напря- 
жения по толщине пленки носит нелинейный ха- 
рактер, и с увеличением числа ^е максимальное 
значение нормального напряжения, которое до- 
стигается на стенке цилиндра, растет. 



Рис. 5. Распределение нормального напряжения по толщи- 
не пленки в зависимости от числа И/е; 1) 50; 2) 100; 
3) 150; 4) 200 


Поскольку при уменьшении \Ѵе толщина слоя 
стекающей жидкости увеличивается, кривая ра- 


спределения нормального напряжения принимает 
все более пологий характер. 
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Рис. 6. Распределение температуры по толщине пленки в за- 
висимости от числа Ре в сечениях: 2 = 0,14-1) 0,0001; 
2) 5; 3) 10; 4) 20; 5) 25; 2 = 0,75 ~ 6) 0,0001; 7) 20 

Численное исследование поля температур пока- 
зало, что при выбранных значениях чисел Ке=5- ІО 4 , 
Рг=3-10“ 7 и числа Пекле (Ре), указанного на рис. 6, 
число Прандтля принимает действительно большие 
значения. Из представленного рисунка видно, что 
вблизи орошаемой поверхности образуется темпе- 
ратурный слой, в котором происходит основное из- 
менение температурного поля. При этом толщина 
этого слоя занимает примерно третью часть от тол- 
щины стекающего слоя вязкоупругой жидкости. В 
этом слое температура вдоль пленки изменяется го- 
раздо медленнее, чем в ее поперечном направлении, 
и с увеличением числа Ре наблюдается уменьшение 
градиента температуры по толщине пленки. 

3. Выводы 

Получено численное решение задачи о стекании 
пленки нелинейной вязкоупругой жидкости по бо- 
ковой поверхности вертикального цилиндра. Пока- 
зано влияние параметров течения и реологических 
характеристик жидкости на профили скоростей и 
температур, а также установлено их влияние на тол- 
щину стекающей пленки при заданных значениях 
объемного расхода. Установлено, что вблизи оро- 
шаемой поверхности образуется температурный 
пограничный слой, в котором происходит основное 
изменение температурного поля. Толщина этого 
слоя занимает примерно третью часть от толщины 
стекающего слоя вязкоупругой жидкости. 
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